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In the context of translationally-invariant insulators, Wilson loops describe the adiabatic evolution 
of the ground state around a closed circuit in the Brillouin zone. We propose that the Wilson-loop 
eigenspectrum provides a complete characterization of the inversion-symmetric topological insulator. 
Through the Wilson loop, we formulate a criterion for nontriviality that indicates a Z classification of 
ID inversion-symmetric insulators. If the ground-state wavefunctions at momenta and 7r transform 
under different representations of inversion, we find that a subset of the Wilson-loop eigenvalues are 
robustly quantized to -1; we identify the number of -1 eigenvalues as a topological index N irl) G Z. 
Physical interpretations of iV(_i) are provided in holonomy and in the geometric-phase theory of 
polarization. In addition, we identify iV(.i) with the number of protected boundary modes in the 
entanglement spectrum. In 2D, we identify a relative winding number W which provides a Z 
classification of 2D inversion-symmetric insulators. For insulators with nonzero W, their Wilson- 
loop eigenspectra exhibit spectral flow that is protected only by inversion symmetry. Hence, W is the 
inversion-analog of the first Chern class C\ (for charge-conserving insulators) and the Z2 invariant 
S (for time-reversal invariant insulators). Finally, we establish relations between the topological 
invariants (W, Ci, E) and the Wilson-loop eigenvalues at symmetric lines in the 2D Brillouin zone. 

PACS numbers: 74.20.Mn, 74.20. Rp, 74.25. Jb, 74.72.Jb 



There is strong evidence supporting the view that the 
topological properties of a condensed-matter system are 
encoded in its ground state aloneP®' From this point of 
view, a translationally-invariant topological insulator^ is 
characterized by a local gauge redundancy - its ground 
state is invariant under a unitary transformation in the 
subspace of occupied Bloch bands. Since all topological 
quantities must be invariant under this transformation, 
a natural object to investigate is the eigenspectrum of 
a Brillouin-zone Wilson loop (W)P^^In the context of 
topological insulators, W describes the adiabatic, unitary 
evolution of n occ occupied Bloch bands around a closed 
circuit in the BZ. More precisely, W is defined as the 
path-ordered exponential of the Berry- Wi lczek- Zee con- 
nection C^ m (k„) = / d d r ul*{r)d kfi uf(r)^^ 

W = Texp[- JdkuCkiku)], (1) 

where uj?(r) = u%(r + R) is the periodic component 
of the Bloch wave. Here, T denotes path-ordering and 
kn G {fci,...,fcrf} denote momenta components in a 
d-dimensional BZ. As Zak demonstrated for n oac — 1, the 
Abelian W is nothing less than the Berry-phase f actor 
acquired by a Bloch wave along a cyclic evolution! 62 * 63 ! 
For general n occ , ([I]) forms a matrix representation of a 
holonomy, i.e., a parallel transport map; holonomies are 
known to have diverse applications in physicsP^HsU 

The Wilson loop is known to provide a characteriza- 
tion of the time-reversal-breaking Chern insulat or^ a nd 
the time-reversal-symmetric (TRS) Z2 insulator! 57 * 67 ! In 
this paper, we employ the Wilson loop to study topolog- 
ical insulators with the simplest point-group symmetry: 



inversion HSEHHIll] Specifically, we analyze how inver- 
sion (I) symmetry manifests in the KV-eigenspectrum 
of Chern insulators and Z2 insulators. In addition, we 
identify topological invariants that are protected only 
by X symmetry, and offer physical interpretations of 
these invariants in holonomy and in the geometric-phase 
theory of polarization. 

We first present a summary of our findings. For 
the ID Z-symmetric insulator, we derive a mapping 
between the I eigenvalues of the ground state and the 
eigenvalues of W; in particular we show that the only 
allowed W-eigenvalues are +1,-1 and complex-conjugate 
pairs; in Tab. |TJ we show an example of this mapping 
for a generic insulator with four occupied bands. We 
propose that the nontrivial I-symmetric insulator is 
characterized by W-eigenvalues fixed at -1; the number 
(N ( _ 1} ) of degenerate eigenstates at -1 is a Z topological 
index, which has the following physical interpretations: 
(i) if we adiabatically transport an insulating ground 
state with n aaa occupied bands around the ID BZ, 
A ( _ 1) of n occ electrons pick up Berry phases of tt (ii) 
in a periodic ID lattice, N (rl) number of maximally- 
localized Wannier functions localize in the middle of 
each adjacent pair of atomic sites; we call this the 
mid-bond site. Moreover, we shall identify N { _ 1} as the 
number of protected mid-gap modes in the entangle ment 
spectru m 1 5 * 70 ^ of an inversion-symmetric insulator! 48 * 71 * 
The presence of these mid-gap modes ensure that the 
entanglement entropy can never be tuned to zero by 
any gap- and symmetry-preserving transformation - this 
has been proposed as a criterion for nontriviality that 
unifies all known topological insulators.^ While these 



mid-gap modes are boundary states, the Wilson loop is 
a bulk property - together they provide consistent and 
complementary characterizations of the Z-symmetric 
insulator. 
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TABLE I. For an insulator with four occupied bands, we 
tabulate the X eigenvalues of the occupied bands at symmet- 
ric momenta {0, 7r} and the corresponding W-spectrum. We 
collect the X eigenvalues fj at each symmetric momenta into 
(£1 £2 £3 £4). + (— ) refers to an eigenvalue of +1 (-1); A A* and 
fi /i* refer to complex-conjugate-pair W-eigenvalues. 



For the 2D I-symmetric insulator, we represent the 
quantum-mechanical electron density by a lattice of 
classical line charges. These line charges extend in y 
in the manner of a ID Bloch wave, but localize in x 
in the manner of a ID Wannier function (WF), hence 
we call them hybrid WF's. By evaluating Wilson loops 
at constant k y , the phases {${k y )} of the W-spectrum 
represent the cc-coordinates of these hybrid WF's; we 
refer to them as Wannier centersP^HZU Upon applying 
an electric field in y, k y may be thought of as a time- 
coordinate in the adiabatic evolution. {$(k y )} represent 
the real-space trajectories of the Wannier centers in x, as 
a function of time k y . We identify the first Chern class, 
C\ 6 Z, as the net ^-translation of the Wannier centers 
in a period 27rp2HsU the TR invariant, S e Z 2 , is the 
change in time-reversal polariza tion over a ha lf-period; 
S is odd in the nontrivial caae. ™ 52 ' 67 ' 69 ' ™ ^ 

Insulators with nontrivial C\ or S exhibit spectral flow 
in the energy spectrum: on a lattice with boundaries, 
a set of protected edge states interpolates across the 
energy gap, such that the conduction and valence 
bands are spectrally connected.^ For these insulators, 
spectral flow is also manifest in their single-particle 
entanglement spectra, where states close to eigenvalue 
1 (0) form the entanglement analog of conduction 
(valence) bandsP^^ Spectral flow is also known to 
exist in the W-spectrum of these insulators, i.e., 
the Wannier trajectories interpolate across the full 
unit cell: ■& € [0, 27r], where $ is the phase of the 
W-eigenvalues; the symmetries that protect spectral 
flow are charge conservation (for the Chern ins ulator) 
and time-reversal symmetry (for the Z2 insulator) ) 57 l 67 l 86 l 
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FIG. 1. Wannier center flows for four different phases. Each 
figure contains three unit cells in an effective ID-lattice along 
x; we track the real-space trajectories of Wannier centers (blue 
circles) in x as k y is varied through two adiabatic cycles, (a) 
An insulator with Chern number Ci = 2 and zero relative 



winding (W); this is modelled by the Hamiltonian (271 with 
parameters m = 3 and 8 — 2. (b) An insulator with relative 
winding W = 1 and zero Ci; this is modelled by Eq. ( |27[ ) with 
m = 3 and S = 1. (c) W = 2 and C\ = 0; this is modelled by 
Eq. (31 1, (d) W — and Ci = 0; the phase has a nontrivial 
polarization, and is modelled by Eq. 1 24 1 with parameters 
a = -1.5, = 1.5,5 = -I. 



It has been demonstrated in Ref. |48]that spectral flow, 
as protected only by I symmetry, does not exist in the 
physical energy and entanglement spectra; I-symmctric 
topological insulators generically have protected mid- 
gap entanglement modes, without spectral flow. In 
this paper, we report that I-protected spectral flow 
manifests only in the W-spectrum, and is characterized 
by a relative winding number W; W provides a Z clas- 
sification of 2D Z-symmetric insulators. Two Wannier 
trajectories are said to wind relative to each other if they 
(i) intersect the same atomic site at a symmetric time 
K y , then (ii) separate and intersect adjacent mid-bond 
sites half a period later (K y + ir). W is defined as the 
number of stable pairs of relatively-winding trajectories. 
For illustration, we plot the Wannier flows for the cases 
(i) Ci = 2, W = (Fig. [l}a), (ii) W = 1, Ci = (Fig. 
n}b) and (iii) a TRS insulator with W = 2 and 5 = 
(Fig. [l}c). While these three nontrivial examples exhibit 
various forms of spectral flow, the absence of spectral 
flow does not imply a trivial insulator. A case in point 
is the 'weak' X-symmetric topological insulator of Fig. 
[T]-d, which may be thought of as a coupled array of ID 
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Z-symmetric topological insulators; here, spectral flow is 
absent (W = C\ = 0), yet the insulator has a nontrivial 
polarization in x. 

As exemplified in Tab. ||J the mapping between I 
and W-eigenvalues also applies to symmetric lines in 
the 2D BZ, e.g., k y = and tt. For symmetric Wil- 
son loops at constant K y — {0,7r}, let us define the 
number of +1,-1 and complex-conjugate-pair eigenval- 
ues as N (tl) (Ky),N ( _ 1) (K y ) and N (cc) (K y ) respectively; 
these quantities constrain the possible Wannier trajec- 
tories that interpolate between K y = and 7r, and by 
implication they also constrain the possible topological 
invariants: W, C\ and S. These constraints are summa- 
rized in Tab. [TTJ where we define Nd as the maximum of 
two quantities: 

N d = max { N (rl ,(ir) - N ( . t) (0) - iV (cc) (0), 

N^M-N^W-N^fr)}. (2) 





£1(1) 


TRS 


W > 


2W + \d\ = N d 


W= \N d ; 
S = W mod 2 
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\Ci\ > N d ; 

Gx = N ( _ 1} {0) - N^(ir) 
mod 2 


S=f(jV ( . 1) (0)-iV ( . 1) ( 7 r)) 
mod 2 



TABLE II. Relations between relative winding W £ Z, the 
Chern number Ci £ Z, the TR invariant 3 £ Z2, and eigen- 
values of the Wilson loop at symmetric momenta. Columns: 
(7(1) denotes a generic insulator with charge-conservation 
symmetry; TRS denotes a time-reversal symmetric insulator. 
N d is defined in Eq. (2). 

In the remainder of the introduction, we extend our 
discussion of the ID I-symmetric insulator. What does 
it mean for a topological insulator to be nontrivial, and 
what role does inversion symmetry play? We define a 
nontrivial insulator by what it is not, up to gap- and 
symmetry-preserving transformations: an atomic insula- 
tor - a lattice of spatially-isolated atoms. How do we dis- 
tinguish the W-spectrum of a trivial insulator from one 
of a nontrivial insulator? In the atomic limit, all Bloch 
eigenfunctions are independent of crystal momentum, 
therefore the connection C vanishes and W equals the 
identity in the occupied subspace. Having W-eigenvalues 
that deviate from +1 is not a sufficient condition for non- 
triviality - by a gap- and symmetry-preserving transfor- 
mation of the ground state, we may find that W is homo- 
topically equivalent to the identity. A case in point is the 
ID time-reversal symmetric (TRS) insulator, for which 
the W-spectrum consists of Kramers degeneracies which 
may be tuned to any value, including +lJ13ffi3 However, 
suppose a subset of the W-spectrum is fixed by symme- 
try to a special value other than +1 - our inability to 
transform the W to the identity indicates that the insu- 
lator is topologically distinct from the atomic insulator. 



In this paper we demonstrate that the discrete nature of 
the simplest point-group symmetry (I) is reflected in a 
discrete W-spectrum - a subset of I-symmetric insula- 
tors have W-eigenvalues that are fixed at -1, hence these 
insulators satisfy the criterion for nontriviality. In ID, 
inversion I maps the spatial coordinate: x — > -x, where 
the inversion center (x = 0) is chosen to lie on an atomic 
site. In the momentum representation, I maps k — > -k, 
and an I-symmetric insulator is characterized by the X- 
eigenvalues of its occupied bands at symmetric momenta 
and 7r. If the wavefunction at and 7r transform under 
different representations of I, we find that a subset of the 
W-spectrum is robustly quantized to -1. To quantify a 
change in the group representations between and 7r, we 
look to the absolute difference in the number of same- 
symmetry bands between and tt: if there are n ( _j(0) 
negative-X bands at and n ( _)(7r) of them at 7T, we find 
that 

tf^K-^Oj-nc-jMI, (3) 

where is the number of symmetry-protected -1 

W-eigenvalues. In principle, N^ 1} is detectable by 
holonomy: if we adiabatically transport an insulating 
ground state with n occ occupied bands around the ID 
BZ, iVj.j) of n occ electrons pick up Berry phases of n, 
which is detectable by interference. In practice, the 
coherent transport of electrons across the entire BZ is 
impeded by scattering from phonons, impurities and 
many-body effects. Despite these diffic ulties, Bloch 
oscillations and Wannier-Stark ladder states^ ^^ have 
been observed in semiconductor superlattices^ ^ and 
with cold atoms in accelerated optical lattices EShli2l 
These systems are the likeliest candidates to realize such 
an interference experiment. 

While the Wilson-loop index N ( _ ±) characterizes the 
bulk of a ID I-symmetric insulator, it is instructive to 
reformulate N { _ ±) from the perspective of the insulator's 
boundary modes. On a lattice with spatial boundaries, a 
variety of topological insulators may be classifi ed thro ugh 
a stability analysis of their boundary modes ISSliill p or 
example, if one boundary mode is stable but two is not, 
we deduce that the classification is Z2. While ID I- 
symmetric insulators have no protected boundary modes 
in the physical energy spectrum, Turner et al. have 
demonstrated that these insulators manifest protected 
boundary modes in the entanglement spectrum - these 
modes localize at ficti tious boundaries created by a spa- 
tial entanglement cut! 46 * 47 * As demonstrated in Ref. H51 
and [711 the number (x £ Z) of stable modes at each en- 
tanglement boundary is equal to the absolute difference 
in the number of same-symmetry bands between and 
7r: 

X=K-,(0)-n ( _,(7r)|. (4) 

Upon proving the identity (|3j) , we identify \ — -^(-i)i 
and thus establish a correspondence between the bulk 
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of a ID Z-symmetric insulator and its entanglement 
boundary. 



The Hamiltonian for a single particle is 



Let us offer a third perspective through the geometric- 
phase theory of polarization, where the phases of 
the W-eigenvalues represent spatial coordinates of the 
maximally-locali zed W annicr functions (WF), or just 
Wannier centers.^ 14 * 115 ! An eigenvalue of +1 implies the 
presence of a Wannier center at each atomic site in the 
periodic ID lattice; -1 implies the presence of a Wannier 
center at each mid-bond site. In a stu dy o f I-symmetric 
insulators with a single occupied bandp^l Kohn demon- 
strated that there are two classes of WF's: 

(i) The WF is either symmetric or antisymmetric about 
the atomic site [x — 0). This occurs when (i-a) both 
Bloch waves with symmetric momenta k = and 7r are 
nodeless at x = 0, or (i-b) they both possess nodes at 
x = 0. 

(ii) The WF is either symmetric or antisymmetric about 
the mid-bond site (x — a/2). This occurs when only one 
of the two symmetric-momentum Bloch waves has a node 
at x = 0. 

Subsequently, Zak connected these observations to the 
modern theory of polarization, and showed that the 
Abelian W is quantized to +1 in case (i) and to -1 in case 
(ii) ) 62 l 63 l Remarkably, (i) and (ii) are distinguished only 
by the group representations of the two Bloch waves at 
and 7r. One aim of this paper is to extend this Abelian 
result to the non- Abelian W, so that we may characterize 
I-symmetry insulators with multiple occupied bands; 
many realistic materials defy a single-band description, 
due to interband degen eracie s such as those enforced by 
time-reversal symmetry! 48 l 86 l 

The outline of this paper: in Sec. |TJ we introduce 
Brillouin-zone Wilson loops and explain their connec- 
tions with polarization and holonomy. We analyze the 
ID I-symmetric insulator in Sec. [II] Here, we (i) derive 
a mapping between X and W eigenvalues, and (ii) for- 
mulate the topological index N { _ 1} which provides us a Z 
classifi catio n. We analyze the 2D I-symmetric insulator 
in Sec. Ill Here, we (a) describe how the W-spectrum is 



H = 



constrained by I and TRS, (b) identify a relative winding 
number W which characterizes Z-protected spectral flow, 
and (c) elucidate relations between various topological in- 
variants (W, Ci,3) and the W-eigenvalues at symmetric 
lines in the 2D BZ. We conclude in Sec. IIVI 



INTRODUCTION TO BRILLOUIN-ZONE 
WILSON LOOPS 



The Brillouin-zone Wilson loop is pedagogically intro- 
duced; we describe its role in holonomy (Sec. |IA| an d in 
the geometric-phase theory of polarization (Sec. |IBj ). In 
Sec. |I C[ we explain the construction of a coarse-grained 
Wilson loop from a tight-binding Hamiltonian, for the 
purpose of numerically computing topological invariants. 



2m 



V(r). 



(5) 



If V(r) = V(r + R) for any lattice vector R, this Hamil- 
tonian is symmetric under discrete translations. Con- 
sequently, H decouples into representations labelled by 
the crystal momentum fc; an eigenstate in the n'th band 
may be written in Bloch form: ipk( r ) — e%kr u k( r )i where 
u k( r ) — M fe( r + -^) i s a function that is periodic in lattice 
translations, and also satisfies: 



(p + hk) 2 

2m 



V(r) 



uUr) = 0. 



(6) 



Each eigenstate has a corresponding projection 
■p£(r, r') = Ufc(r) v%*(r'); the many-body ground- 
state is a single Slater determinant of all single-particle 
eigenstates with energies less than the Fermi energy. 
The topological properties of an insulator are invari- 
ant under transformations of the Hamiltonian that 
preserve both the energy gap and the symmetry that 
stabilizes the topological phase. We perform one such 
transformation by setting the energies of all eigenstates 
m£ below (above) the Fermi energy equal to e_ (e+). 
Denoting the projection onto the n occ occupied bands 
as V%~° = Y^i=i ^fei we express the resultant flat-band 
Hamiltonian as 



H F (k) = (e. 



-)n 



■e + I. 



(7) 



Eq. (fn) has a gauge redundancy which is not appar- 
ent in wok - the ground-state projection V^" is invariant 
under alocal U(n 0{:c ) gauge transformation in the Tri- 
dimensional subspace of occupied bands: 



with m, n = 1 . 



and A^ 1 =Ml 



A. Wilson Loop as arising from Holonomy 

The adiabatic transport of a ground state at initial 
momentum to a final momentum k involves a uni- 
tary rotation of the basis vectors u 1 ^ in the subspace of 
occupied bands. This U(n occ ) rotation is affected by a 
Wilson-line matrix W fc ^ fc (i) that maps the subspace of 
occupied bands at fc (s) to the subspace of occupied bands 
at k. W is known to satisfy a parallel transport equation 



with the connection C defined as 



(8) 



CT(*p) = J d d rut{ry^ut{r). (9) 

Here, k^ S {/ci, . . . , k^} denote momenta components in 
a d-dimensional BZ. Eq. $h is pedagogically derived in 
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App. [X] This differential equation has the path-ordered 
solution 



k«-k(0 



(10) 



for a path C that connects momenta fc and k {i) . (10) 



is a matrix representation of a holonomy, i.e., a parallel 
transport map!^ If k = fc (l) modulo a reciprocal vector, 
£ forms a closed circuit in the BZ and we denote the 
resultant U(n occ ) Wilson loop as W; c/. (JlJ. 



B. Wilson Loops in the Geometric-Phase Theory 
of Polarization 



Let us derive a well-known relation between W and 
the polarization of a ID insulator. For a group of n oaa 
occupied bands, the delocalized Bloch waves ip™ form an 
orthonormal basis in the occupied Hilbert space; m = 
1,2,..., n oce . Assuming a periodic gauge (ip™ = '0fc+27r)' 
we formulate the theory of polarization in an alternative 
basis of localized Wannier functions (WF) : 



n occ 

ikR 

n=l 



o(k) jn r k {x). (ii) 



Each WF *Cj)(a; - #) i s labelled by the unit cell ReZ 
and an index j = 1,2, ... n occ . If n occ = 1, O(k) is 
a momentum-dependent phase; if n occ > 1, 0(k) is a 
U(n occ ) matrix that affects rotations in the subspace of 
degenerate bands. The gauge freedom in O(k) is fixed, 
up to trivial U(l) phase windi ngs, by requiring that the 
WF's are maximally localized!*^ Equivalently, we re- 
quire that the WF's are eigenfunctions of the projected 
position operator V° cc xV occ ^^ 



-pooc £ -poc 



27T 



R) ¥ j) (x~R) =0. (12) 



Here, V "" projects to the occupied subspace, and all 
lengths are in units of the lattice constant. We show 
in App. [c] that the spectrum coincide with the 

phases of the W-spectrum. Through (12 1 we map the 



electron density to classical point charges at the posi- 
tions {??w)/27r + -R}; we call these positions the Wannier 
centers. If d = 0, there exists a Wannier center at each 
atomic site; if = n, there exists a Wannier center in 
the middle of two adjacent atomic sites - we call this the 
mid-bond site. For an atomic insulator, W equals the 
identity in the occupied subspace; this implies that all 
the Wannier centers sit at the atomic sites. 



C. The Tight-Binding Wilson Loop 



orbitals which diagonalize the atomic Hamiltonian. The 
Hamiltonian of ^ reduces to the variational form 



H = ^2 C L i h ( k )} a fs c kfj 



(13) 



with orbital indices a, /3; [h(k)] a g are the matri x eleme nts 
of (p5| ) in the Bloch basis of Lowdin orbitals P^t^ j n 
( |13| and the rest of the paper, we sum over repeated 
indices. Let us denote the j-th normal mode operator as 



7 fc 



[Ujl]p c\p\ in bra-ket notation, the corresponding 
projection is V k = \U 3 k ){U 3 k 



the tight-binding connection A as 



7j fc |0>(0| 7j - fe . We define 



(14) 



The tight-binding Wilson loop W A is defined as the path- 
ordered exponential of the tight-binding connection A: 



W A 



Texp[- J 



(15) 



A discretized expression of W A is obtained by divid- 
ing a BZ loop C into infinitesimally-separated momenta: 
{k m +G, fc (JV) , k {N - 1} . . . fc <2) , fc (1) , fc <0) } with N > 1. Defin- 
ing Vj,° c = Y^=i T^l as tne projection to the occupied 
bands, W A may be expressed as a path-ordered prod- 
uct of projections, sandwiched by tight-binding eigen- 
functions at the base and end points: 



N 



[w A {c)r n = (u% 0)+G \ run: |c^ 0) >. 



(16) 



The product of projections are path-ordered along £, 
with the earlier-time momenta positioned to the right. 
In App. [D] we show that A is related to the continuum 
connection C, as defined in through 



Ckik^) = A k (kfj) - i [f M ]k + iR ll 



(17) 



where f M is the position operator in the basis of occu- 
pied bands and R^ is the average atomic displacement 
in the direction [i. Since A and C differ by [f^]k — R^, 



we expect that the tight-binding Wilson loop (15) and 
the continuum Wilson loop ([T]) generically have differ- 
ent eigenspectra. However, we show in App. [E] that 
both Wilson loops are identically constrained by inver- 
sion symmetry, and consequently their spectra are nearly 
identical. In particular, the topological index A r ( _ 1) , as de- 
fined in the Introduction, may be extracted from either 
Wilson loop. To simplify the presentation in the next 
Section, we speak only of the tight-binding Wilson loop, 
which we henceforth denote as W . 



As defined in ([!]) , the continuum Wilson loop W is ex- 
pensive to compute numerically. In this Section, we in- 
troduce the tight-binding Wilson loop W A , which is com- 
putable with relative ease. In the tight-binding approx- 
imation, we restrict our attention to nearly-degenerate 



II. WILSON-LOOP CHARACTERIZATION OF 
THE ID INVERSION-SYMMETRIC INSULATOR 

We highlight distinctive features of the I-symmetric 
Wilson loop (Sec. II A), and present a mapping between 
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the W- eigenvalues and the X eigenvalues of the ground 
state (Sec. II B I. In Sec. II C we formulate the topologi- 



cal index iV ( _ 1} £ Z that classifies Z-symmetric insulators. 
In addition, we relate N ( _ 1} to a well-known Z2 index that 
distinguishes the electric responses of these insulators. 



The symmetry constraints (19 1,(20 1 and (21 1 are derived 



in App. [EJ these results are applied in the next Section, 
where we derive a mapping between W-eigenvalues and 
the I eigenvalues of the occupied bands. 



A. Constraints on the Wilson Loop due to 
Inversion Symmetry 

As shown in App. [Dj a tight-binding Hamiltonian with 
X symmetry satisfies ph(k)p = h(-k), where p is the 
overlap matrix: 



Pa/3 



dx 



(18) 



a,/3 are orbital indices and {4> a } are eigenfunctions of 
the atomic Hamiltonian. This implies that the occupied 
eigenstates t/™ at ±fc are related through inversion by a 
U(n oce ) 'sewing matrix' B^: 



\U r _, 



mn* 

k Pa/3 



US 



kid- 



(19) 



n=l 



Here, m, n are indices that label the occupied bands. 
Due to the U(n acc ) gauge freedom in the occupied 
subspace, the sewing matrix B^ is gauge- variant. In 
deriving the gauge-invariant VF-eigenspectrum, we will 
employ a gauge in which t/2„v are eigenstates of X 
at inversion-invariant momenta k l " v ; B k i nv is then a 
diagonal matrix with diagonal elements equal to the X 
eigenvalues of ?72 nv . 

Given B k and its analytic properties, we derive the 
following properties for an X-symmetric Wilson loop: 

(i) W is unitarily equivalent to its Hermitian adjoint, i.e., 
the set of W-eigenvalues is equal to its complex conju- 
gate: 



{expii?} = {exp-zt?}. 



(20) 



This implies that the eigenvalues of W are constrained 
to ±1 or otherwise form complex-conjugate pairs. 
Equivalently stated, the Wannier centers localize (i) at 
the atomic site (+1) or (ii) at the mid-bond site (-1) or 
(iii) form pairs that localize equidistantly on opposite 
sides of each atomic site (AA*). In all three cases, the 
spatial configuration of Wannier centers is invariant 
under a spatial inversion x — > -x. 

(ii) Let £ be a BZ loop beginning at fc (0) = -it and end- 
ing at it. The Wilson loop along C may be decomposed 
into two Wilson lines that each connect two symmetric 
momenta: W(£) = Wn-oWo^-i. Up to a change in 

orientation, W^-^o is mapped to Wo< n by an inversion 

k — > -k. It follows from the sewing condition (19) that 



W(£) = [Wo^l'BoWo* 



(21) 



B. ID: Mapping between Wilson-loop and 
Inversion Eigenvalues 

For a subset of I-symmetric ground states, their W- 
spectra comprise only ±1 eigenvalues. This occurs if the 
I eigenvalues of occupied bands at either k = or it are 
the same; then the W-spectrum comprises the diagonal 
elements of the product of sewing matrices: Bq B^ . 
Proof: If at k — 0, the occupied bands have X eigenvalue 
£, i.e., the sewing matrix Bq — £/, then from (21) we 
deduce that W = B B^. If instead B^ = £ /, then W is 
equivalent to Bq B n by a unitary transformation. 

Let us define a quantity n s that equals zero for the 
above subset of ground states; n s > if the occupied 
bands have nonidentical X eigenvalues at both and 7r. 
Definition: For the occupied bands of an insu- 
lating phase, let us define the number of positive 
and negative X eigenvalues, at fc mv = {0,7r}, by 
n (+) (fc lnv ) and n ( _)(fc mv ). Given this set of four num- 
bers {n( + )(0),n(_)(0),n( + )(7r),n(_)(7r) }, we identify the 
smallest of the four and label it as n s , i.e., n s is the 
number of the Fewest Bands of One Symmetry among 
both symmetric momentum fc ,nv ; we call these bands the 
FBOS. We label the momentum where FBOS lies as k s 
and the X eigenvalue of FBOS as £ s . Let us identify the 
FBOS for the following examples. 

(a) Consider a two-band insulator with X eigenvalues 
(++) at k = and (+-) at k = it. The FBOS are 
the negative-I bands at k s — 0. None exists, so n s = 0. 

(b) Suppose we had a four-band insulator with X eigen- 
values (+ H ) at k = and (+ + H — ) at k = it, the 

FBOS is the single negative-X band at k s = it, so £ s = -1 
and n s = 1. 

(c) If a subset of the four numbers {n (±) (fc inv )} are 
equally small, we may denote any number in this subset 
as n s . In a two-band example, we may encounter 
rc. (+) (0) = ?i(-)(0) = n {+) {it) = n^^it) — 1. Then, one 
may label any of the four possibilities as the FBOS. 

Mapping: Given an inversion-symmetric insulator that 
is characterized by the quantities { n (±) (fc' nv ), n s , k s , £ s }, 
its Wilson-loop eigenspectrum consists of: 

(i) (n (+) (fc s + it) — n s ) number of -£ s eigenvalues, 

(ii) (n ( _ } (k s + it) — n s ) number of £ s eigenvalues, and 

(iii) n s pairs of complex-conjugate eigenvalues. 

In the above examples, the W-spectrum of insulator (a) 
comprises one +1 and one -1 eigenvalue; for insulator 
(b), there are one +1 eigenvalue, one -1 eigenvalue, 
and one complex-conjugate pair; insulator (c) has one 
complex-conjugate pair only. The proof of this mapping 
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is detailed in App. [Gj The interested reader also may 
refer to App. [FJ where we undertake the case studies 
of the one- and two-band W's; these case studies offer 
an intuitive understanding of the above mapping, and 
also provide an alternate derivation that is specific to 
one and two occupied bands. For an insulator with one, 
two and four occupied bands, we tabulate the possible 
mappings in Tab. |IIf[ |IV| and |T| respectively. 

Let us define a Z index, N { _ 1} , as the number of -f 
eigenvalues in the W-spectrum. It is possible that one 
or more pairs of complex-conjugate eigenvalues are ac- 
cidentally degenerate at -1; we exclude them from the 
definition of N ( _ 1} . Through the above mapping, we de- 
duce that N ( _ 1} is the absolute difference in the number of 
same-symmetry bands between and n, and quantifies 
the change in the group representation; cf. Eq. In 
the following Section, we argue that a nonzero is an 
indication of topological nontriviality. 



X eigenvalues 


W-spectrum 


{(+) (+)} 


[+] 


{(+) HI 


H 



TABLE III. For an insulator with one occupied band, we tab- 
ulate the X eigenvalues of the occupied band at symmetric 
momenta {0,7r} and the corresponding W-spectrum. + (— ) 
refers to an eigenvalue of +1 (-1). {(+) (— )} may refer to ei- 
ther (i) a positive-I band at k = 0, with a negative-Z band 
at 7r, or (ii) a negative-Z band at 0, with a positive-! band at 
7r. If two sets of X eigenvalues (from two distinct insulators) 
are related by a global change in sign, they are mapped to 
the same W eigenvalue. For example, both {(+) (+)} and 
{(— ) (— )} are mapped to W = +1. 





X eigenvalues 


W-spectrum 


(i) 


{(++) (++)} 


[++] 


(ii) 


{(++) (+-)} 


[+"] 


(iii) 


{(++) (— )} 


[— ] 


(iv) 


{(+-) (+")} 


[AA*] 



TABLE IV. For an insulator with two occupied bands, we 
tabulate the X eigenvalues of the occupied bands at symmet- 
ric momenta {0, n} and the corresponding W-spectrum. We 
collect the X eigenvalues at each symmetric momenta into 

(6 6). 



C. ID: Topological Invariants from the Wilson 
Loop 

1. Z Index: Number of Robust -1 Wilson-Loop Eigenvalues 

A nonzero N { _ 1} index obstructs W from being tuned 
to the identity - the atomic limit; this is a sufficient 



condition for nontriviality. While complex-conjugate W- 
eigenvalues [A, A*] also deviate from the atomic value of 
+1, they are not a sufficient condition for nontriviality. 
The exact value of A is not hxed by symmetry; this ar- 
bitrariness reflects the range in this equivalence class, 
i.e., it is possible to tune the Hamiltonian and sweep the 
interval of allowed A while preserving both gap and sym- 
metry. In particular, we can always tune the complex- 
conjugate eigenvalues to the trivial limit of [+1,+1]. It 
is possible that pairs of complex-conjugate eigenvalues 
are accidentally degenerate at -1. Unlike the N {rl) val- 
ues of -1, these extra degeneracies are not protected by 
T - they generically destabilize under a soft deformation 
of the ground-state, hence they do not indicate a non- 
trivial phase. Hence, symmetry dictates that is the 
minimum number of -1 eigenvalues; there are two impli- 
cations: 

(i) Suppose we begin with a nontrivial insulator with 
N (rl) > 0, and we would like to transform it to an 
atomic insulator while preserving I symmetry. The en- 
ergy gap must close a minimum of -/V M) times at a sym- 
metric momentum, before the atomic limit = 0) is 
reached. In each gap-closing event, a pair of opposite-sign 
I-eigenvalues are inverted between occupied and unoccu- 
pied bands, thus reducing by one. 

(ii) As formulated in the entanglement spectrum, the 
criterion for nontriviality is the existence of protected 
boundary modes in the gap of the single-particle en- 
tanglement spectrum. For an X-symmetric insulator, 
the number (x € ^) °f stable modes at each entangle- 
ment boundary is equal to the absolute difference in the 
nu mber of same-symmetry bands between and tt; cf. 
Q 1481711 Thus, we identify x with the W-index N ( _ 1)7 and 
show that the two formulations of nontriviality are equiv- 
alent. 



2. Z2 Polarization Index: Determinant of W 

The determinant of W, which we dehne as T>, is the 
exponentiated polarization of the ID insulator. Since all 
W-eigenvalues are either ±1 or form complex-conjugate 
pairs, T> is quantized to ±1 - the classification of the 
electric response is Z2, as is recognized in works such as 
Ref. US] and [S3] Moreover, T> is only determined by the 
number of -1 eigenvalues: T> = (-l)^ 1 '. Let us relate T> 
to the I eigenvalues of the ground state. From ^ we 
have that 

n occ 

V = (_i)K-)(o)-«c-)«l = Yl ]] 4™., (22) 

fe inv =0,7r m=l 

where £J?? nv is the I eigenvalue of the m'th band at sym- 
metric momentum k lm . This concludes our discussion for 
ID. 
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FIG. 2. The Wannier center flows for a 2D insulator with X 
symmetry. Each figure contains three unit cells in an effective 
ID-lattice along x; the positions of the atomic sites are indi- 
cated by black dots. The ground state contains two occupied 
bands, hence there are two Wannier centers in each unit cell. 



(a) This phase is realized in the model Hamiltonian ( 24 1 , with 
parameters a = -1.5, ft = 1.5 and 8—1. Upon k y — > k y + 2tt, 
the pair of Wannier centers in each unit cell exchange posi- 
tions. There is no net transfer of charge between unit cells, 
hence Ci = 0. (b) This phase is realized in the same model 
with a — -1.5 and f5 = 0. Upon k y — > k y + 2ir, one Wannier 
center adiabatically flows to its neighboring unit cell on the 
left. This represents a quantized Hall current, hence C\ = -1. 



III. WILSON-LOOP CHARACTERIZATION OF 
THE 2D INVERSION-SYMMETRIC INSULATOR 

The Wilson loop W is known to encode the first Chern 
class G\\ we present a summary of this relation in Sec. 
|III A| In Sec. |III B[ we impose I symmetry and inves- 
tigate how the symmetry constrains W and the allowed 
Chern numbers. W is further constrained if the insula- 
tor is also time-reversal symmetric (TRS) - this is ex- 
plored in Sec. |III C| In Sec. |IIID| we introduce a rel- 
ative winding number W that characterizes insulators 
with I-protected spectral flow. In addition, we derive 
relations between W, C\, the time- reversal invariant S, 
and Wilson-loop eigenvalues at symmetric lines in the 2D 
BZ. 



A. Wilson Loops and the First Chern Class 

A well-known relation exists between Wilson loops 
and the integer quantum Hall effect. Let us apply an 
electric field along y, which adiabatically translates all 
single-particle states in the parameter space k y . To 
probe for a quantum Hall response, we study the k 
evolution of the ground-state polarization (in x) 
interested in Brillouin-zone W's at constant k y : W ky — 



v 

we are 



Texp(- J dk x Ck(k x )). Let us denote the n occ eigenval- 
ues of W fe by the set {exp )}. The geometric phase 



i?™ /2-7T represents the center of a hybrid Wannier func- 
tion (WF), which extends in y in the manner of a ID 
Bloch wave, but localizes in x as a ID WF; we refer to 
these phases as the Wannier centers.^! The derivative of 
the geometric phase, $™ = d$ m /dk y , is interpreted as 
the real-space velocity (in x) of the m'th Wannier center 
at time k y . By integrating the velocities of all Wannier 
centers over a period 27r, we obtain the net quantum 
Hall current. Thus, we identify the first Chern class as 
the center-of-mass winding: 



(23) 



As we have shown in Sec. |IB| polarization is directly 
related to the continuum Wilson loop, which is defined in 
Q. However, the winding number in the tight-binding 
Wilson loop, as defined in Sec. |I C| is identical to that in 
the continuum Wilson loop. This follows because their 
connections differ by an operator that is periodic in k y 
(cf. Eq. (17) ). For the purpose of computing Chern 



numbers, both Wilson loops give identical results. 
For illustration, we consider the 4-band model 

h(k) = i r 13 + f (cos k x + cos k y ) (r 30 + r 03 ) 



(ri 



31 



2 

sin k. 



(r 2 



32, 



sin k„ — T, 



0.3 



(24) 



+ | (cos k x ) (cos k y - 5) (r 03 - r 30 ), 



where h(k) is a matrix in the tight-binding basis; cf. ( |13[ ). 
Tij are defined as ai <g) Tf, ao (to) is the identity in spin 
(orbital) space; 0^=1.2,3 (^=1,2,3) are Pauli matrices in 
spin (orbital) space. The Hamiltonian possesses an I 
symmetry: To 3 ft-(fc) To 3 = h(-k). The Fermi energy is 
chosen so that there are two occupied bands in the ground 
state. We tabulate the 1 and W-eigenvalues for various 
choices of the parameters (a, j5 and 8) in Tab. [V] In 
Fig. [2j-a we plot the W ky -spectrum for a trivial insulator 
(a = -1.5,/? = 1.5,(5 = 1); Fig. [2]-b corresponds to a 
nontrivial insulator with C\ — -1 (a — -1.5, (3 = 0,5 = 0). 



Parameters 


X eigenvalues 


Wa' h eigenvalues 


a 


P 


5 


(0,0) 


(JT.0) 


(0,7T) 


(n,n) 


Ky=Q 




-1.5 








(++) 


(++) 


(++) 


(+") 


[++] 


[+-} 


-1.5 


1.5 


1 


(++) 


(++) 


(+-) 


(+-) 


[++] 


[AA*] 


-1.5 


1.5 


-1 


(++) 


(+") 


(++) 


(+") 


[+-} 


[+-] 



TABLE V. The X and Wk v eigenvalues of the ground state 
of Hamiltonian (|24|), for various parameters. 



B. 



The Inversion-Symmetric Wilson Loop and the 
Integer Quantum Hall Effect 



Let us investigate the spectrum of the Z-symmetric 
W k . As derived in App. ITTJ we find that W k is equiva- 
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lent to the Hermitian adjoint of W_ k by a unitary trans- 
formation, i.e., the sets of W-eigenvalues at ±k y are 
equal up to complex conjugation: 



(a) 



{exp id ky } = {exp -id. ky } , (25) 
At K y £ {0,tt}, the ID 



as may be verified in Fig. 2 
line of states behaves like a ID X-symmetric insulator in 
two respects: (i) the eigenvalues of Wk are constrained 
to ±1 or otherwise form complex-conjugate pairs, (ii) 
In ID, the X eigenvalues at k = and it are related 
to W-eigenvalues through the mapping of Sec 



MB 



2D, the X eigenvalues at momenta (0,K y ) and (n,K y ) 
are related to the eigenvalues of Wk v through the same 
mapping. 

Let us define the number of robust -1 eigenvalues in 
the spectra of Wo and W^ as N ( _ i:i (0) and A r ( _ 1) (-7r) re- 
spectively. During the adiabatic evolution, N^^Ky) is 
the number of Wannier centers that localize at each mid- 
bond site at time K y . A difference in the indices N { _ 1} (0) 
and iV ( . 1 )(7r) implies a net Hall current; moreover, the 
parity of the Chern number is determined through 

^(-1,(0) - AVdM = Ci mod 2. (26) 

The two parities of C\ correspond to the following situ- 
ations: 

(i) Suppose the parities of N { _ 1} (K y ) differ. Between 
Ky = and 7r, an odd number of Wannier centers must 
interpolate between the mid-bond sites (W-eigenvalue of 
-1) and the non-mid-bond sites, which include (a) the 
atomic sites (+1) (see Fig. [2]-b), and (b) the complex- 
conjugate sites (A A*). The net translation of Wannier 
centers in the interval k y £ [0,7r] is half an odd integer. 



It follows from (23) and (251 that C\ is odd 



(ii) An analogous argument emerges when the parities of 
N(_ r) (K y ) are equal. Now an even number of Wannier 
centers must interpolate between the mid-bond sites (at 
time K y = 0) and the non-mid-bond sites (at time K y = 



The 



7r). One possible scenario is illustrated in Fig. 
conclusion is that C\ is even. 

It follows from Eq. p6]), ^ and ([22} that the product 
of all X eigenvalues (over all occupied bands at every 
symmetric momenta) has the same parity as C\. 



C. The Inversion- and Time-Reversal-Symmetric 
Wilson Loop 

Our aim is to highlight distinctive features of Wilson 
loops with both X and time-reversal symmetries; we illus- 
trate these features with a four-band model Hamiltonian: 

h(k) — (2 — m — cos k x — cos k y ) To3 + 5 sin k y T\2 

+ sinfc,, (r 31 + rn) + smky (r 2 i + r 02 ), (27) 

with matrices IYy defined as Oi <g) Tf, a (r ) is the 
identity in spin (orbital) space; 0^=1.2,3 (^=1,2,3) are 



ky/7l 




- 1 a/Tc 1 



FIG. 3. The Wannier flow of a 2D insulator with both X and 
TR symmetries, (a) This phase has a trivial TR-invariant; 
it is realized in the model Hamiltonian (271, with parameters 



m = 4.3, S — 0. (b) This phase has a nontrivial TR-invariant; 
it is realized in the same model, with m — 3, 8 = 0. 



Pauli matrices in spin (orbital) space. The Hamiltonian 
is I-symmetric: To3 h(k) To3 = h(-k). The Fermi energy 
is chosen so that there are two occupied bands in the 
ground state. We tabulate the X and W-eigenvalues for 
different choices of parameters m and S in Tab. |VI[ In 
this Section we set 5 — 0, so the Hamiltonian is also 
time-reversal symmetric (TRS): T/i(fc)T _1 = h(-k), with 
T = iT2o K; K is the complex-conjugation operator. 
The two classes of TRS insulators are distinguished by 
a Z2 invariant S, which is the change in time-reversal 
polarization over half an adiabatic cycle; S is odd for 
the nontrivial classP^ In Fig. (|3]-a) and ([3]-b), we have 
plotted the W fc -spectra for both Z2-trivial (m = 4.3) 
and nontrivial (m = 3) phases respectively. 



Parameter 


T eigenvalues 


Wk v eigenvalues 


m 


5 


(0,0) 


(tt,0) 


(0,7T) 


(7T,7r) 


Ky=0 


Ky = TV 


3 





(++) 


(++) 


(++) 


( — ) 


[++} 


[ — ] 


3 


1 


(++) 


(++) 


( + + ) 


(— ) 


[++] 


[ — ] 


3 


2 


(++) 


(++) 


( + + ) 


(— ) 


[++] 


[ — ] 


4.3 





(++) 


(++) 


( + + ) 


(++) 


[++] 





TABLE VI. For various choices of the parameters m, 5 in the 



Hamiltonian (27 1, we write the corresponding (a) X eigenval- 
ues at the four symmetric momenta and (b) the eigenvalues 
ofWjr,. 

As derived in App. |TJ TRS imposes the following con- 
straints on the W fc spectra: 

(i) The sets of eigenvalues at ±k y are equal, i.e., 

{expi$ ky } = {expitf_ fcy }. (28) 

(ii) The Ws at symmetric momenta satisfy 

wL = e- 1 w Ky e, (29) 
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FIG. 4. (a) An Z-symmetric phase with relative winding 
W = 2; spectral flow is protected by I symmetry alone. This 
is the ground state of Hamiltonian (31 1, (b) Spectral flow is 



interrupted with an I-breaking perturbation. 



with an antiunitary operator that squares to -I. This 
implies that every eigenstate of VVk v has a degenerate 
Kramer's partner, (i) and (ii) imply that if one Wannier 
center produces a Hall current Ih, its Kramer's partner 
produces a time-reversed current that cancels Ih- As 
shown in Rcf. [57] and 1671 the Z 2 invariant may be 
extracted from Fig. ^ in the following manner: in the 
region k y £ [0, 7r], d £ [-7T, 7r], let us draw a constant-?? 
reference line at any value of d. If the Wannier trajecto- 
ries intersect this reference line an odd number of times, 
the phase is nontrivial, and vice versa. 

By imposing T symmetry as well, we arrive at the fol- 
lowing conclusions: 

(a) Due to I symmetry, the Wk -spectra at K y = {0, 7r} 
consist of ±1 and complex-conjugate pairs; the addi- 
tional constraint of Kramer's degeneracy implies that the 
spectra is composed of pairs of [+1,+1], pairs of [-1,-1] 
and complex-conjugate quartets [A A A* A*]. Since time- 
reversal and I commute, the two states in a Kramer's 
doublet must transform in the same representation un- 
der I - this limits the possible X eigenvalues in a TRS 
ground state. 

(b) From p5l) and (p8| we derive 



{expi$ ky } = {exp-itf_ fey } = {expit?_ fcy }, 



(30) 



which indicates that the flow of the Wannier centers in 
one quadrant, say $ £ [0, 7r] and k. 
the flow in the full range, £ 



v 

7T, 7T 



£ [0, it], determines 



and k y £ 



D. Inversion-Protected Spectral Flow and the 
Relative Winding Number 

The Chern insulator and the TRS topological insulator 
exhibit spectral flow in the W-spectrum; the symmetries 
that protect spectral flow are, respectively, charge conser- 
vation and time-reversal symmetry. In this Section, we 
report spectral flow of a third kind, which is protected 
by T symmetry alone. In our first example, we consider 
the eight-band model 

h{k) = (-1-cos k x -cos k y ) 6 3o + sin k x (Q 3W + ©no) 



sinky (6 2 io + ©c 



0.8 sin k x (6311 +8111), (31) 



with matrices Qijk defined as <Ji ® Tj <8> 7fc; 0^=1,2,3 are 
Pauli matrices in spin space; for i,j — {1,2,3}, <£> jj 
are products of Pauli matrices in a four-dimensional 
orbital space; ctq (tq ® 70) is the identity in spin 
(orbital) space. This Hamiltonian is I-symmetric: 
0030 h(k) 0030 = h(-k), and time-reversal symmetric: 
T/i(fc)T- 1 = h(-k); here T = i@ 2 ooK, and K imple- 
ments complex conjugation. The Fermi energy is chosen 
so that there are four occupied bands in the ground 
state. We tabulate the I and W-eigenvalues in Tab. 
VII and also plot the W ky -spectrum in Fig. Hj-a. With 
TRS, C\ — 0. The change in time-reversal polarization 
over half an adiabatic cycle is 2 (even), hence the 
TR invariant is trivial™ Yet, Wannier trajectories 
interpolate across the full unit cell: i9 £ [0, 27r). Let 
us softly break I symmetry, while maintaining TRS, 
with the perturbation: 0.4 cosk x Q022 + 0.4 cosk y ©112- 
As evidence that spectral flow persists only with I 
symmetry, we find in Fig. [4]-b that the spectrum is now 
gapped. 



X eigenvalues 


y^K v eigenvalues 


(0,0) 


M) 


(0,7T) 


(n,n) 


Ky=0 


Ky=1T 


+ + ++ 


+ + ++ 


+ + ++ 




+ + + + 





by reflections. This is illustrated in Fig. ([3]), where each 
quadrant is bounded by dotted lines. 



TABLE VII. For the Hamiltonian (|3TJ, we write the (a) 2 
eigenvalues at the four symmetric momenta and (b) the eigen- 
values of Wj{„ . 



I-protected spectral flow is characterized by a nonzero 
relative winding number W, which is defined in the fol- 
lowing way. Two Wannier trajectories are said to wind 
relative to each other if they (i) intersect the same atomic 
site at a symmetric time K y , then (ii) separate and inter- 
sect adjacent mid-bond sites half a period later (K y +tt). 
W is defined as the number of stable pairs of relatively- 
winding trajectories. To make this definition precise, we 
outline a procedure to identify W: 

(a) Count the number of Wannier trajectories that di- 
rectly connect mid-bond and atomic sites in the quad- 
rant {d £ [Q,ir],k y £ [0, 7r]}; call this number n\. By 
a direct connection, we mean a smooth trajectory that 
flows without interruption. We consider three examples: 
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in Fig. [5]-a (blue) and[5]-b (blue), n\ = 0; m = 2 in Fig. 
Ha. 

(b) Count the number of trajectories that directly con- 
nect mid-bond and atomic sites in another quadrant 
{d £ [—it, 0),k y £ [0, 7r]}; call this ri2- In all three exam- 
ples, n\ — n 2 - 

(c) W is the minimum of {ni,^}. We find that the 
insulator of Fig. |3[a has relative winding W = 2; W = 
in the other two cases. 

For each pair of relatively-winding trajectories, one of 
the pair ($i(k y )) has winding number +1 on the torus 
{i9 G [-7r,7r),fcj, <E [-7T, 7r)}, and the other has winding 
-1. In principle, it is possible that l Qx{k y ) has winding 
2n+l (n € Z + ), while its partner $2 (fey) has winding 
-2n-l. However, n > implies that the trajectories 
dxiky) and $2 (fey) cross at a non-symmetric momentum 
(fey + 1 {0,7r}); such degeneracies are not protected 
by symmetry - by a ground-state deformation that 
preserves both the energy gap and T symmetry, we may 
turn crossings into anti-crossings and thus reduce n to 
0. Since W is the number of stable relatively-winding 
pairs, we eliminate all such accidental degeneracies 
before carrying out the above procedure to identify W. 
All accidental degeneracies fall into two categories: (i) 
at non-symmetric momenta, there may be accidental 
crossings of two or more trajectories, as illustrated in 
Fig. [5|b (red). A slight deformation results in level 
repulsion, and turns crossings into anticrossings (Fig. 
[5jb (blue)). (ii) At symmetric momenta {0,7r}, we 
rule out complex-conjugate-pair eigenvalues that are 
degenerate at either the atomic or mid-bond site. In the 
example of Fig. [5]-a (red) , there is one such degeneracy 
at the atomic site when K y = 0, and another at the 
mid-bond site when K y — 7r; upon perturbing the 
Hamiltonian, this degeneracy splits, as shown in Fig. 
[§a (blue). 

Let us derive a sufficient condition for relative wind- 
ing, and simultaneously relate W, C\ and the eigenvalues 
of W at constant K y — {0, n}. We define N (tl) (K y ), 
N(-i)(K y ) and N (cc) (K y ) as the number of +1,-1 and 
complcx-conjugate-pair eigenvalues of W K , respectively; 
N {tl) {K y ) + N irl) (K y ) + N (ac) (K y ) = n oec ,\he number of 
occupied bands. Define the quantity Nd as in Eq. ([2]); if 
Nd > 0, there are Nd trajectories that directly connect 
the mid-bond site (at K y ) and the atomic site (at K y +ir); 
otherwise, there are none. Of these Nd trajectories, one 
or more pairs may relatively wind, and the rest wind with 
center-of-mass motion, thus contributing to the Chern 
number C\. The parity of G\ is constrained as in (26 1. A 



sufficient condition for relative winding is that \C\ \< Nd, 
in which case 2W = N d - \C X \. If \d\ > N d , then 
W = 0. For 1- and time-reversal-symmetric (I+TRS) 
insulators, C\ = 0, hence the relative winding is related 
to Nd through: 2W — Nd- These relations, as summa- 
rized in Tab. [TTJ also imply that W is isotropic, in the 
following sense. We define {(p(k x )} as Wannier trajecto- 
ries of the Wilson loop at constant k x ; if {^(fcy)} exhibits 
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FIG. 5. Figures in red provide two examples of accidental 
degeneracies. Upon a gap- and symmetry-preserving pertur- 
bation, these accidental degeneracies split, as shown in blue. 



relative winding W, then so will {ip(k x )}. This claim is 
substantiated in App. [j] 



1. Relative Winding of Insulators with both Inversion and 
Time-Reversal Symmetries 

In this Section we study the relative winding of insu- 
lators with both I and TRS (I+TRS); we shall relate 
the relative winding W with the TR invariant S. While 
both W and 3 characterize Wannier trajectories with 
no center-of-mass motion, they differ in many important 
respects. For I+TRS insulators with nonzero relative 
winding, the parity of W determines the TR invariant: 
S = W mod 2; 2 is odd in the nontrivial class. To prove 
this, we apply the rule: modulo 2, S equals the number 
of Wannier trajectories that intersect a constant-t? 
reference line.^ With I+TRS, only one quadrant, e.g. 
{$ € [0,n], k y € [0,7r]}, is independent. In the rest of 
this section, we denote coordinates in this quadrant 
by ($, k y ). Two cases are possible: (i) W number of 
trajectories directly connect points ($, k v ) = (tt, 0) and 
(0, 7r), or (ii) W trajectories connect (7r,7r) and (0,0). 
If n occ = 2W, there are exactly W intersections with 
the reference line, hence S = W mod 2. If n occ > 2W, 
it is possible in case (i) that: (i-a) an extra trajectory 
connects (0,0) and (0,7r), (i-b) a trajectory connects 
(ir, 0) and (tt, n), (i-c) if there exists a complex-conjugate 
quartet [AiAiA^A^] in the spectrum of Wrt =oj a pair 
of trajectories may connect (0, n) with the complex- 
conjugate site at (Ai,0), and (i-d) if there exists a 
complex-conjugate quartet [A2A2A2A2] in the spectrum 
of Wif =7r , a pair of trajectories may connect (tt, 0) with 



12 



the complex-conjugate site (A2, tt). In all scenarios, these 
extra trajectories intersect an even number of times 
with the reference line - the parity of the number of in- 
tersections is decided by W alone. The proof is complete. 

While only the parity of W matters to the Z2 classifi- 
cation under TRS, W provides a Z classification under T 
symmetry, and hence a more complete characterizaton. 
This distinction may be understood from a stability 
analysis of the Wannier centers. Since Kramer's degener- 
acy is two-fold, four or more Wannier centers generically 
experience level repulsion. Consider for example the 
W = 2, I+TRS model of Fig. [4}a. Sitting at the atomic 
site (at K y — n) are four Wannier centers which are 
constrained by Z-symmetry - they do not experience 
level repulsion. If we now break X-symmetry while 
preserving TRS, these four Wannier centers destabilize 
and split to form two pairs of Kramer's doublets, thus 
breaking spectral flow; see Fig. [4}b. 

For I+TRS insulators with zero relative winding, it is 
possible that spectral flow is completely absent and the 
insulator is trivial. However, an I+TRS insulator with 
four or more occupied bands may have a nontrivial TR 
invariant without relative winding. To distinguish these 
two cases, we look to the indices iV ( . 1} (0) and N^^tt). 
Since W — 0, all N^^Ky) Wannier centers that origi- 
nate from the mid-bond site (at K y ) must flow to either 
a complex-conjugate site or a mid-bond site (at K y + tt). 
By conservation of trajectories, A/ M) (K y ) — iV ( _ 1} (K y + tt) 
number of trajectories must connect the mid-bond site 
(at K y ) to complex-conjugate sites (at K y + tt); here we 
have assumed N ( _ 1) (K V ) > N ( _ l} {Ky + tt). Since C\ = 0, 



we know from (|26|) that Nr.^Kyj 
within one quadrant, e.g., 



even 



- N ( _ 1} (K y + tt) is 

{4 £ [0,7T], ky G 

[Ky ,K y +n}}, (JV ( . 15 (K y )-N w (Ky+n) )/2 number of tra- 
jectories connect the mid-bond site (at K y ) to complex- 
conjugate sites (at K y + tt). Two cases arise: (i) if 
(N(.i)(K y ) — N^(K y + 7r))/2 is odd, there must be at 
least one other trajectory that interpolates between K y 
and K y + tt. Then the sum of all trajectories that con- 
nect to complex-conjugate sites is even, as required by 
Kramer's degeneracy. Hence, TRS enforces a zig-zag pat- 
tern of Wannier flows, as illustrated schematically in Fig 
[6}a and -b. (ii) For even (N^ (K y ) ~ N W {K V + tt))/2, 
Kramer's degeneracy is satisfied without additional tra- 
jectories, and the resultant Wannier flows are gapped, as 
in Fig. [6}c and -d. Therefore, 



S = § (AVdW-AVdM) mod 2 
mod 2; 



= 3^ 



(32) 



for the last equality, we applied the definition of AT d in 
Eq. |2j) a nd that N {ca) (K y ) is a multiple of four; c.f. 
Sec. Ill C The relation ( 32 ) applies to I+TRS insulators 



with W > as well, since we have proven 2 = W mod 
2 in this Section, and previously identified 2W = iV^ 
in Sec. HID The relations between S, W and Nd are 




FIG. 6. Schematic illustrations of wannier flows with zero rel- 
ative winding. In figures (a) and (b), |(-/V ( _i)(0) — (tt)) = 
1 (odd); the Wannier trajectories follow a zig-zag pattern, (c) 
|(A+i)(0) — A+i)(7r)) = (even), and the Wannier flows are 
gapped, (d) §(iV(-i)(0) - A^dW) = 2 (even). 



summarized in Tab. HH 

Consider for example the eight-band model: 



h(k) = ( -2 - cos k x - cos k y ) <d 033 + | 
+ 5 sin fc^ (9ii 3 + 6313 
+ ( 6 sin k x + I sin k y ) ( 9 



^023 



g cos k y 0ooi 

e 2 i ) 

e 2 i 3 ) + e 3o, (33) 



§ sin k y (9 20 



with matrices 9yfc defined in (|3 1|). This Hamiltonian has 
the same symmetries as that in ( |3l[ ), namely I and TRS, 
and the ground state is defined to be the four lowest- 
energy bands. We tabulate the I and W-eigenvalues 

[7}a, 



in Tab. VIII 



As illustrated in Fig. [TJ-a, the relative 
winding is trivial, but partner-switching occurs with help 
from the complex-conjugate quartet at K y = 0. In this 
W = model, spectral flow is protected by TRS alone. 
As evidence, we deform the Hamiltonian with a TRS- 
breaking term, Ah(k) = 8 sin k x ( 9 22 + 0120 + 9m ) + 
2 sin k y 9 123, which preserves both the energy gap and T 
symmetry. As illustrated in Fig. [7]-b, Kramer's degener- 
acy is now lifted - the quartet [AAA* A*] at K y — splits 
into two separate doublets [Ai A*J, [A2 A?>]. This contrasts 
with a previous example, where we broke TRS in a 3- 
nontrivial phase; as shown in Fig. [l]-b, the resultant W- 
spectrum is not gapped, due to a nonzero relative wind- 
ing. 



X eigenvalues 


Wk v eigenvalues 


(0,0) 


M) 


(0,7T) 


(n,n) 


Ky=0 


Ky=TT 


+ H 


+ + — 


+ H 




A A A* A* 


+ H 



TABLE VIII. For the Hamiltonian (|33|, we write the (a) X 
eigenvalues at the four symmetric momenta and (b) the eigen- 
values of Wjf„ . 
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FIG. 7. (a) An Z-symmetric phase with zero relative wind- 
ing; spectral flow persists due to TRS alone. This is the 
ground state of Hamiltonian (33 1, (b) Spectral flow is inter- 
rupted with a TRS-breaking perturbation. 



Appendix A: Derivation of parallel transport 
equation ([8]) 

Let us define it£ as eigenfunctions of the flat-band 
Hamiltonian "H F (fc), with energy e_; the Hamiltonian is 
derived in Eq. Q. n is a band index that runs over 
the occupied bands, and k labels the crystal momentum. 
u k form basis vectors in the subspace of occupied bands. 
Upon adiabatic evolution through the Brillouin zone, the 
state uj (j ) at some initial momentum fc <l) (at time tt) is 
mapped to a different state v k at final momentum k (at 
time t). In this process, the state acquires a dynamical 
phase J. e.dt. In addition, there is a unitary rotation of 
basis vectors in the occupied subspace; this rotation is 
affected by a U(n occ ) Wilson-line matrix W^ki')'- 



«2(r) 



-if* S.dt 



E «rwwrA« (ai) 



IV. CONCLUSION 



In the adiabatic approximation, v k remains an eigenstate 
of the T-L F (k) with energy e_. Let us apply the time- 
dependent Schrodinger equation and sum over repeated 
indices: 

e. v n k (r) = H t (k) v n k (r) = id t [e"< 4 F - dt u™(r) H&>] . 

(A2) 



In this paper we have developed a topological classifi- 
cation of inversion-symmetric insulators (ISI) based on 
the Wilson loop. We propose to classify the ID ISI with 
an index iV ( _ i; £ Z, which coincides with the number 
of prote cted boundary modes in the entanglement 
spectrum! 48 ! 71 ! In 2D, we identify an invariant W that 
characterizes spectral flow in the Wilson-loop spectrum; 
W provides a Z classification of the 2D ISI. In contrast, 
W has no analogs in both the entanglement and physical 
energy spectra, where there exists no spectral flow that 
is protected by inversion symmetry alone. A theory 
to formulate N { _ 1} and W as topological obstructions 
to certain gauges remains to be developed; the experi- 
mental verification of these invariants also remains an 
open question. Inversion is but the simplest point-group 
symmetry - we may ask how the Wilson-loop spectrum 
is influenced by other point-group symmetries, e.g., C n . 
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Upon cancellation of the dynamical phase and replacing 
dt = Ylu=i (dk^/dt) d kfl for a c?-dimensional BZ, we ar- 
rive at Eq. pi). 



Appendix B: Properties of Wilson loops 

In this Appendix we derive several properties of Wilson 
loops, which apply to both continuum and tight-binding 
versions, as defined in Sec. |I A| and |l C| respectively The 
derivations have been carried out with the tight-binding 
connection A k (cf. Eq. (14) ) and for a ID BZ; they 



are trivially generalizable to the continuum connection 
Cfe (cf. Eq. Q ) and to higher dimensions, 
(i) Let us denote a Wilson line, over a path with start 
point k (i) and end point fc (/) , as <-&(«• The Wil- 

son line between two infinitesimally-separated momenta 
is W^_"^ fe = (Uj£_ e \U%). Let us define loop C as tracing 
a path between base point k m and end point k m + 2tt; 
we divide the C into infinitesimally-separated momenta: 
{fc (0) + 2tt, k (N) ,k (N - 1} . . . fc <2) , fc (1) , fc <0) }, with N » 1. Let 
^fc = Em=i |^™)(^fci denote the projection to the oc- 
cupied bands at momentum k. The Wilson loop may be 
discretized: [W(C)] mn 

- w ml w lq w op w pn 



TV 



=(% +2I |Tn^ K(o>). 



(Bl) 
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The product of projections are path-ordered (symbolized 
by T) along C, with the earlier-time momenta positioned 
to the right. 

(ii) Let C be a path that connects two arbitrary momenta 
k {1> and k {2> , and C T be the same path with opposite path- 
ordering. The Wilson line satisfies the unitary condition 



w fe ( 2) ^ (1 ,(e)t = w, 



(B2) 



For k {2> = k lr> + 2tt, the above relation generalizes to 
Wilson loops. From (B2) one can derive that for a fixed 



loop, the eigenspectrum of W is independent of the base 
point. 

(iv) Under a U(n occ ) gauge transformation, the prod- 
uct of projections in Eq. (Bl) is invariant, but the 



eigenf unctions at the base and end points transform as 
fe) -> E m \ u ?(o)) M mn (*<°>). Hence, the Wilson 
loop transforms as W -> M^(k m ) W M(fc <0) ), and the 
eigenspectrum of W is gauge-invariant. 



Appendix C: Wilson Loops and the Projected 
Position Operator 

Let us demonstrate that the phases of the W-spectrum 
coincide with the eigens pect rum of the projected posi- 
tion operator V oco x V co EE1 Here, P occ projects to the 
occupied subspace of the translationally-invariant Hamil- 
tonian (15]), which have Bloch eigenf unctions ip k (x) = 
e ikx u n ~yy e are interested in eigenfunctions "3/ that 
satisfy 



{V° a °xV aa °-±$) *(x)=0 



(CI) 



for some eigenvalue $/2tt. We expand if? in the subspace 
of occupied Bloch waves: 



".CI r. 

*(*) = E / i fnk 

n=l J 



(C2) 



In the periodic gauge (-0J? = "0^+2^)' tne ac tion of the 
projected position operator on the wavefunction / may 
be decomposed into an in tra-band operator d k and an 
inter-band operator CP^I 



(M\v°~xv°°°\*)=i^+i x: P(fc)r m / mfe , 

m— 1 

(C3) 



where C is the non-Abelian connection that is defined in 
Q. In the general solution to (IClj), the wavefunctions 
at two different momenta (k and are related by the 
Wilson line 



-/,*(*) C(g)dg-|mn 



(C4) 



It follows from the periodic boundary condition on / that 

n occ 

E Ink = f m ,k+2-n = ^ f m . k , (C5) 

hence we have shown that the eigenspectrum of 
■pace £ -pocc coineideg w ith the phases of the eigenspectrum 
of VV. Furthermore, let us derive the eigenfunctions of 
•pocc £-pocc Tj enne the U(n occ ) matrix Q(k) such that its 
columns are the eigenstates of the Wilson loop at base 
point k: W fe+27r <_fc. That is, W k+2 ^ k = Q{k) DQ^(k), 
where D is a diagonal matrix that contains the eigenval- 
ues of W. While the eigenvalues of W do not depend on 
the base point k, the eigenfunctions Q(k) do. The matrix 
Q(k) is related to the matrix Q(k' 7^ fc) by a Wilson line: 

=> Q(fc') = wL fc 'Q( fc ) = Wfe'^Q( fc )- (C6) 

We label the j'th diagonal element in D by exp ), 
where j runs from 1,2 ... rt occ . If i?«)/2tt is an eigen- 
value of 'P occ & •p occ j so is any integer addition to /2tt, 
as is consistent with i?vJJ being a phase. For each occu- 
pied band, there exists an infinite ladder of eigenvalues: 
$ { j> ) /27r = i?W/2tt + G Z. The index i? labels the 

unit cell where the eigenfunction is localized. We will 
choose the convention that -n < fid) < 71-. Let us then 
label the eigenfunctions of "P occ x "P occ by j = 1 . . . n occ 
and £ Z. The wavefunction & oc Q(k) n j so as 

to satisfy the eigenvalue equation (C5|. In addition, we 



multiply the j'th column of the eigenmatrix Q{k) by a 
momentum-dependent phase, so as to ensure periodicity: 

[f^U.k+2, = e-^+^n/^ Q( k + 2n) nj 

= e -«(^*r)*gV*r J- W ^ Q( fc)m . 
m— 1 

= e -»(*+^ 5 /*r e »w q (% . = [f U) ]nk _ (C7) 
In the second equality, we have applied the relation |C6| . 



In summary, to each eigenvalue 1?^ /2ir there corresponds 
an eigenfunction of 'P oco x V° cc : 



^(x-R) = ^ 



n=l 



27T 



Q(fc) ni ^(x). (C8) 



Appendix D: The Tight-Binding Formalism and the 
Relation between the Full and Tight-Binding 
Connections 

In the tight-binding variational approximation, the 
Hilbert space is reduced to n tot atomic orbitals (p a (r — R) 
on each lattice site (R) ; <\> a are eigenstates of an atomic 
Hamiltonian i? atom , and a = 1,2, . . . ,n tot is an or- 
bital/spin index. Due to the finite overlap of atomic 
orbitals on different lattice sites, we are motivated to 
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construct an orthonormal basis that preserves the point- 
group symmetrie s of <j> a . Such a basis is realized with 
Lowdin functions 111ZH11SI which are defined as 



0,0 



with the hermitian overlap integral 



(Dl) 



Kf = I dd r ~ Ri) M x ~ R i) = A T*- (° 2 ) 



/3a* 



Here, i,j are indices for lattice sites; we sum over re- 
peated indices. The orthogonality of Lowdin function 
reads as 



J d d rip* a (x-R)cp (x-R')=d a0 S R , Rl . (D3) 



Let us define the operator I by its action on functions: 
I f(r) — f{-r). If the atomic Hamiltonian preserves I 
symmetry, i.e. [H atom ,I] = 0, then one may define a 
Hermitian, unitary, n tot x n tot sewing matrix p in the 
basis of atomic orbitals: 



P a p= I d d r<j>* a (r)I<f> fj (r) = (<j> a \i\<p ). (D4) 



The Lowdin functions ip, which are constructed from <f> 
through (Dl), transform identically under point-group 
symmetry operations liiZJ hence I <p a (r) — <p a {-r) — 



<Pg(r) p Sa as well. In Hamiltonians with discrete transla- 
tional symmetry, we form the Bloch sums 



R 



e- ik ^ip a (r - R), (D5) 



which are periodic in lattice translations. The tight- 
binding Hamiltonian (cf. ( 13 ) ) is defined as 



[h(k)] a0 = (u ka \ H(k) \u kfj ) 

d«r<p* a (r-R)[£ + V(r)]<p (r) (D6) 



In the tight-binding approximation, the periodic compo- 
nent of the Bloch wave ip% is 



1 



Ra 



^r-R) {U n ]a(pa(r _ R)) (D7) 



where [U£] a is the unitary eigenmatrix that diagonalizes 
(D6). Let us relate the full connection Ck (cf. |9])) to 
the tight-binding connection Ah (cf. (14)). By inserting 
(D7) into (pj), we obtain (17), where we have defined ry, 
as the position operator in the basis of occupied bands: 



VJfc 



e iHR-R') 

R.R 1 



x / d d r <p* p (r-B!)r M (p a (r-R), (D8) 

and Rfj, — J2r R^/N is the average atomic displacement 
in the direction fi. One can prove that the connections 
A and C are anti-Hcrmitian, while the position opera- 
tor r M is Hermitian. Under a U(n oca ) gauge transforma- 
tion, U™ —> UlM^ m (k), the matrix [r^]k undergoes a 
unitary transformation [r^k — > M (k)[r ^]kM (k)' 1 - thus, 
the eigenspectrum of is gauge-invariant. 



Appendix E: Analytic Properties of the Sewing Matrix and Constraints on the Wilson Loop due to Inversion 

Symmetry 



Let us employ the bra-ket notation that is introduced in Sec. IC As defined in Eq. (19), the sewing matrix links 
occupied eigenstates U k of the tight-binding Hamiltonian at ±k: B k nn = [t/™]* p a/3 [U%]g = (U™\ p \U%)\ p is an 
overlap matrix that is defined in Eq. ( 18 ). Here, m, n are indices that label the occupied bands. We choose a gauge in 
which, at inversion-invariant momenta, B k ™ is a diagonal matrix with diagonal elements equal to the T eigenvalues 
of J/Snv We prove that Bk is unitary. Let us label the unoccupied bands by the primed index m! = n occ +l, . . . , n tot . 
If the ground state is insulating, i.e., ([/_™ | p \U k ) = for all n occupied bands and all m' unoccupied bands, 



(irz\p\uj;) (u%\ P \uiy 



"tot 



(U l k \ p\U m k ) {U m k \ p\UZ) = (U l k \ p 2 \U%) = 5 ln . (El) 



In the second equality, we have used the identity: (U™ | p \U£} 
relation. Furthermore, the Hermicity of the inversion operation ■ 



0; the third equality required the completeness 
[BkV 



implies B. 



Bk 1 - 



Applying ( 19 ) and 



/, we derive the following condition for a Wilson line between two infinitesimally-separated momenta fc (1) and 
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fc (2) : 



•\Ajmn / TT m \TT n \ jymo Ijto I 2 \ \ r/rd* r/mo r-tt; lol r>nl* 

W-km^-km — \ u -k( 2 ) \ u -km) — D ki 2 ) \ u km I P \ u kw) n kW — D k&) r^w 2 ) o k{1) 



(E2) 



This relation is generalizable to a Wilson line between arbitrary momenta. We divide the finite path between fc (1) and 
k m into infinitesimally-separated momenta: {fc (1) , fc (1) + A, k m + 2A, . . . , fc (2) — A, fc (2) }. Following Eq. (Bl I, 



yVkW^kW - KV fc(2>^fc(2)_A ,/V fc< 2 )-A<~-fc( 2 )-2A ••• KV fe(i)+A^fe(i) 



n a W_ k (i)_ A ^_ k (i) -o.fcfi) 



— D -kV) KV -fc( 2 ><— fcC 2 )+A ^-fcW+A ^-fct^+A KV -fc( 2 >+A^-fc( 2 >+2A i3 -fe( 2 )+2A ■ • • XJ — A 

— £>_ fe (2) » / V_ fe (2)^_ fc (2) +A ' / > / _A;( 2 )+A^-fe< 2 )+2A ' • ■ KV -fc(i)-A^-fe<i) - -°-fc( 2 ) "-feW-f— fcW ^-fcW- 

Applying this result to a Wilson loop between -it and 7r, 



W mn — R mr y\; rs Ejns* _ umr VAifs nns* _ nmr VAlfrs 

' K 7T^ 7r -7T -7T-*— 7T 7T ^ -ir ' v tt4 it *r -<rr ' ' V J— 



77< 7T 7T 



\B 



snlf 



(E3) 



(E4) 



Here we have used the identities (B2) and B k i m - = £?L„. JE4| informs us that W is equivalent to its Hermitian adjoint 



through a unitary transformation, hence we derive (20) for VV. Let us apply (E3) to derive (21) 



In 



(E5) 



We have derived (20) and (21) for the tight-binding Wilson loop ( cf. (15) ), for which the connection involves 
tight-binding eigenfunctions UJ}. Similar results apply to the continuum Wilson loop ( cf. ), for which the 
connection involves u k - eigenfunctions of the Bloch Hamiltonian ([6]). We may similarly define a U(n oac ) sewing 
matrix: B r k nn = J d d r u™* (r) u k (-r) that links Bloch eigenfunctions at ±fc. Within the tight-binding approximation, 
we may relate u k and U k through (D7), and identify B as identical to the sewing matrix B of (19). This identity 
implies that the mapping between I and W eigenvalues, as presented in Sec. |IIB| applies to both tight-binding and 
continuum Wilson loops. If the ground state is characterized by n s = (n s is defined in Sec. II B I, the spectra of both 
Wilson loops are identical, and comprise only ±1 eigenvalues. If n s > 0, the spectra of both Wilson loops comprise 
the same numbers of ±1 and complex-conjugate-pair eigenvalues; they may differ only in the phases of the complex 
eigenvalues. 



Appendix F: Case Studies of the ID 
Inversion-Symmetric Insulator 

1. Case Study: 1 Occupied Band 

With one occupied band, there is only one Wannier 
center per unit cell. In each unit cell, the Wannier 
center can only be at a atomic site or at a mid-bond 
site; only these spatial configurations in a periodic lat- 
tice are invariant under inversion. If the Bloch waves 
at k — and n transform under different represen- 
tations of I, W = -1. Proof: We shall employ no- 
tation that is defined in (21). Being unitary, Wq^-tt 
must be of the form expjvd). (21) informs us that 
W(£) = B-K exp {-id) B exp = B„ B Q , which is a 
product of the 1 eigenvalues of the sole occupied band at 
k = and ir. The proof is complete. This result can be 
verified by a model tight-binding Hamiltonian 



h(k) = — (a + cosfc)r3+ sinfcr2, 



(Fl) 



with Ti defined as Pauli matrices in orbital space. This 
Hamiltonian has an I symmetry: T3 h(k) T3 = h(-k); the 
insulator is trivial when a > 1, and nontrivial when -1 < 
a < 1. 



2. Case Study: 2 Occupied Bands 

For I-symmetric insulators with two occupied bands, 
W is a 2 x 2 matrix; our case study captures many qual- 
itative features of larger-dimensional W's. We assume a 
general form for Wo< -k that satisfies unitarity: 



W 0M = e" 




W = 1- (F2) 



Inserting (F2) into (21), we arrive at W(£) 




c*dti{e -e Q ) 
e(ici 2 £ 2 + i«) 



(F3) 



where ^ inv and £^ inv are diagonal elements of the sewing 
matrix i? fc mv. We exhaust the possible 1 eigenvalues 
{((.mv}, solve the characteristic equations and derive the 
W spectra. Our results are tabulated in Tab. |IV| 

Cases (i)-(iii) of Tab. IV may be summarized as: if 



the 1 eigenvalues of occupied bands at either k = or tt 
are identical, i.e., if either sewing matrix Bq or B^ is pro- 
portional to the identity, then the W-spectrum comprises 
only ±1 eigenvalues; its eigenvalues are the diagonal ele- 
ments of the product Bq B^. 
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Model parameters 


X eigenvalues 


w 

spectrum 


a 


P 





7T 








(++) 


( + + ) 


[++] 


1.5 





(+-) 


( + + ) 


[+-] 


1.5 


-1.5 


(— ) 


( + + ) 


[— ] 


1.5 


1.5 


(+") 


( + ") 


[AA*] 



TABLE IX. A model of an I-symmetric ground state with 
two occupied bands: for various choices of the parameters 
a,/3 in the Hamiltonian ( |F4[ ), we write the corresponding I 
and W-eigenvalues. 



In case (iv) of Tab. IV 



we encounter (a) occupied 
bands with nonidentical I eigenvalues at both k = and 
tt, and (b) a complex-conjugate pair of W-eigenvalues 
(AA*): a pair of Wannier centers are positioned equidis- 
tantly on opposite sides of each atomic site. The ex- 
act position is not determined by symmetry; this arbi- 
trariness reflects the range in this equivalence class of 
I-symmetric insulators, i.e., it is possible to tune the 
Hamiltonian and sweep the interval of allowed A while 
preserving both the insulating gap and I symmetry. This 
implies that case (iv) is trivial. Why? Let us tune the 
Hamiltonian to a limit in which bands of identical repre- 
sentation are fully coupled, i.e., the positive-I (negative- 
I) band at k = -it adiabatically evolves into the positivc- 

I (negative-I) at 0, through the Wilson line Wo< «■ For 

example, if Co = = = _ £o, the two W-eigenvalues 



± 2i|c||d|. Tuning |c| -> and |d| ->■ 1 



are \d\ — 

effectively decouples the two-band W into two one-band 
W's; each Abelian W connects bands of the same repre- 
sentation between k — and tt, hence each contributes 
+1 to the spectrum, and W — > I. As discussed in the 
Introduction, if W is tunable to the identity, then the 
insulator is in the same equivalence class as the atomic 
insulator, hence we conclude that case (iv) is trivial. 

In contrast with complex-conjugate-pair W- 
eigenvalues, protected -1 eigenvalues obstruct W from 
being tuned to the identity. The nontrivial insulators 
are cases (ii) and (iii) of Tab. IV which have one and 
two -1 eigenvalues respectively. We distinguish between 
case (iii) and the degenerate limit of case (iv), in which 
W — > -I; this is the limit in which bands of identical 
representation are fully decoupled, i.e., the positive-I 
(negative-I) band at k = -tt adiabatically evolves into 
the negative-I (positive-I) at 0. In case (iii), the equality 
W = -I is robust against gap- and symmetry-preserving 
transformations of the ground-state, while this is not 
true in case (iv). 

Let us verify our results in Tab. |IV| with the model, 
tight-binding Hamiltonian 

h(k) = -r 03 + o.i r 13 + (r 21 + r 3 i) sin k 
+ \ ( a (r 30 + r 03 ) + p (r 30 - r 03 ) ) cos k, (F4) 

with matrices defined as cr, ® Tf, (tq) is the iden- 
tity in spin (orbital) space; erj =1) 2 l3 (ri=i 2,3) are Pauli 



matrices in spin (orbital) space. The Hamiltonian is I- 
symmetric: ro 3 h(k) To 3 = h(-k). The Fermi energy is 
chosen so that there are two occupied bands in the ground 
state. We tabulate the W spectra for various choices of 
the parameters a and /3 in Tab. IX 



Appendix G: Proof of mapping between I and W 
eigenvalues 

We employ notation that has been defined in Sec. |II B| 

l (±)( fcinv )> FBOS, fc s ,£ s ,? 



As defined in Eq. (19 I, 



B k is a U(n oac ) sewing matrix linking bands at ±fe 
through I. At symmetric momenta fc inv , a gauge is cho- 
sen in which B k ^ is a diagonal matrix with elements 
{(ifi;™), &(k inv ), ■ ■ ■ , S,n occ (k inv )} equal to the I eigenval- 
ues at fc mv . At momentum k s , we pick a convention that 
the fewest bands of one symmetry (FBOS) are indexed 
by to — 1, 2, . . . , n s and the rest of the bands at k s are 
indexed by to = n s +l, . . . , n occ . We evaluate the Wilson 
loop W over a path that begins at k s -TT, sweeps the in- 
terval [fc s -7r, k s + tt] and ends at k s + tt; the W-spectrum 
is independent of the base point, as shown in App. [B] In 
the rest of the section, we simplify notation and assume 
that W = yVk t +w<r-k t -ir- Let us reproduce a result pre- 
sented in (El]): W = B ks+7I B kg Z, with Z defined as 
the Wilson line Wk s ^k 3 -n'- Z^ 1 = Z\ Let us define Y as 
the n nrr x n nrr matrix 



Since Bl +ks 



K= !(! + & B n+ks W). 
/, the matrix elements of Y are 



(Gl) 



(G2) 



z=i 



Here we have applied the unitarity condition Z^ Z = I to 
express [Z^B k Z\ l3 = Yl=l Uh) Z*, Z aj = -&(<Jy - 



2 J2b=i Z bi Zbj)- We deduce from (|G2] that Y is a rank- 
rig projection matrix. If n s — 0, Y is the zero matrix. We 
define V QlQ2 - Qm as m x m submatrices in Y that lie on 
the intersections of rows numbered by {ai, 0.2, ■ ■ ■ , ct m } 
with columns numbered by {ai, a.2, ■ ■ ■ , ct rn }. For exam- 
Pie, Y 1 =E"=i^,i Za,u and 



y2 3 _ / 2^o=l ^a,2 ^a,2 l^ a =l A a,2 ^a,3 

I n s ry* ry \^^.s ry* ry 

\2^a=l Z a, 3 Z a,2 Z^ a =l Z a,3 Z a,3 . 



(G3) 



The determinant ofamxm submatrix Y is also called the 
m x m minor of Y . According to a well-known theorem 
in linear algebra, the rank of a matrix is equal to the 
largest integer r for which a nonzero r x r minor e xists , 
therefore det Y aia2 -- a ™ = if to > n s . Applying (|G1~|), 



the characteristic equation det [A I — W] = is equivalent 
to 



= det[ -is B ka+7T ] det[ XI -W 
= dct[(-£ s \B ks+1T -I) + 2Y] 



(G4) 
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We claim that the determinant in the second line of ( G4 1 is equal to a polynomial in A of order 2n s , multiplied by 
the factor (-£, A-l)[' l <+)( fe s+^)-n s ] ^ X-l)l n (-)( k >+*)-™'] : where U(+) (fc s + tt) + tt)) is the number of positive-X 

(negative-X) bands at k s + tt. Upon proving this claim, we deduce that there are (n (+) (fc s + tt) — n s ) number of -£ s 
W-eigenvalues and (rt ( _)(£; s + tt) — n s ) number of +£ s W-eigenvalues. Furthermore, we apply a result presented in 
(20 1: the W-eigenvalues can either be ±1 or form complex-conjugate pairs A, A*. Thus, the zeros of the polynomial 



R(X) correspond to n s complex-conjugate pairs of eigenvalues. 
Proof : 

We define A; = -£ s d(k s + tt) A. By a binomial-like expansion, we express the characteristic equation as: 



o=E 2 '" E 



det Y aia2 - a " 



n 



(G5) 



{a>l 



The sum ^2{ ai a2 a } runs over the n occ choose m combinations of the indices {ai . . . a m }. An example of the 
characteristic equation for an insulator with 4 occupied bands is 



= (A 1 -1)(A 2 -1)(A 3 -1)(A 4 -1) + 



2F 1 (A 2 - 1)(A 3 - 1)(A 4 - 1) + 8 det [Y } (A 4 -l) + comb. 1 



+ 



4 det [F 12 ] (A 3 - 1)(A 4 - 1) + comb. 2 



16 det [Y 



12:i4i 



(G6) 



where [/(1234) + comb. 1] = /(1234) + /(4123) + /(3412) + /(2341) and [/(1234) + comb. 2] = /(1234) + /(1324) + 
/(1423) + /(2314) + /(2413) + /(3412). Applying the rank-minor theorem, only the first n s +l terms in ( |G5| ) are 
nonzero: 



0=^2 m det F" 1 " 2 -"™ J] (Arl). 

m—0 {Q 1 ,Q 2 ...Q m } i=X\i-^a.x- > OL2—0tm 



(G7) 



Let us organize this expansion. We first consider a term in (G7) with a particular combination of m band indices 



given by the set {ai . . . a m } in the superscript of Y . Each band index ai has a corresponding X eigenvalue £ ai 
at k = k s + tt; in this set we may define m + (m,— ) as the number of positive (negative) I eigenvalues in the set 
{£ai(^s + 7r) . . .£ Qm (fc s + 7r)}; note m + + m_ = m. The presence of each band in this set implies that a factor 
of (-£ s £ai(k s + tt) A-l) is absent in the product Ili^iWan a 2 a (Aj-1)- ^ n * ne remainder of this proof, we use the 
convention that a, (ft) are positive-X (negative-I) band indices at rr+k s . We may organize the expansion by collecting 
terms with the same m , : 



E 2 ™ E 

rn—O m-4- =0 



-m+] 



(6A-1 



5, 



Mi i . m — rn I 



with S x . y defined as the sum of all minors det Y ai ■■■ Qm +^ 1 —0m- with m + — x and m_ = y: 



c - 

— 



{ai-Uj /3i.../3 y } 



By definition, there are n (+) (fc s 



S{ ai ... Qa 



(n ( _)(fc s + 7r)) number of positive-X (negative-X) bands at k s 



(G8) 



(G9) 



hence the 



0i-M 



runs over n (+) (fc s + 7r) choose x combinations of positive-X bands, multiplied by n ( _)(fc s + tt) 



choose y combinations of negative-X bands. As the form of (G8 1 suggests, nonzero terms in the expansion have values 



of m + ranging from a minimum of to a maximum of n s . That the expansion (G7) includes nonzero terms with 
and nonzero terms with m , 



m 



n s , m- 



0, m_ 



is a consequence of our construction: by definition of 



FBOS and n s , both n (+) (k s + tt) and n^{k s + tt) > n s . We observe that: 



(i) Terms in the expansion with m 4 



. m,- — are proportional to (-£ s A-l) 



because n, factors of 



(-£ s A £, ai (k s + tt)-1) with £ Qi — +1 are removed from the product ri ( A^-l) . All other terms in the expansion have 
greater powers of (-£ s A-l). 

(ii) Similarly, terms in the expansion with m + = 0,m_ = n s are proportional to (£ s A-l)[ n (-^ fes+7r ^ _n "] and all other 
terms in the expansion have greater powers of (£ s A-l). 

(i) and (ii) imply that the common factor of all terms in the expansion is (-&, A-l) K+)( fc a +*■)-« J (£, A-l)K-)( fe »+T)-"»]. 
The characteristic equation is thus expressible as 



= (-&A-1 



[n {+) (k s +n)-n s ] 



(6A-1) 



(_)(*i+7r)-n s ] 



i?(A) 



(G10) 
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with R(X) a polynomial of order 2n s - the claim is proven. 



Appendix H: Proof of Eq. (25 | 



Let C ky be a path in the constant-fcj, interval 
[(-7r, k y ), (77, k y )\; we choose the path-ordering convention 
that k x increases along the path C hy . We define (i) the 
time-reversed path T C hy , which sweeps the same inter- 
val with opposite path-ordering (i.e., decreasing k x ), and 
(ii) the I-mapped path £ C ky , which sweeps the interval 
[(-7T, -k y ), (n,-ky)] with decreasing k x . In combination, 
£TC ky = C~ ky \ this observation leads us to the following 
identity: 



W* (£"*») = W{TC ky ) = B^ ky W(£TC ky )B! iT ky 



B^ ky W{C ky )Bl^ 



(HI) 



where B k are sewing matrices defined in Sec. |II A| The 
first equality follows from ( B2 ) , the second from ( E3 1 . 
Since B 7Tjky — B.^^ ky is unitary, we have shown that VV 
along C ky is equivalent to the Hermitian adjoint of VV 
along C~ ky by a unitary transformation, as advertised. 



Appendix I: Wilson Loop of the 2D Time-Reversal 
Symmetric Insulator 

The time-reversal operator is written as T — Q K, 
with Q^ 1 — and K the complex-conjugation oper- 
ator. On spin-half single-particle states, T 2 = — /. which 
implies Q T — — Q. We define a matrix that sews occupied 
Bloch bands at ±k through T: Vjf = (U l _ k \T\U 3 k ); i,j = 
l...n occ . Q T — ~Q implies V k T — -V. k . By a similar 
proof as one presented in App. |EJ one can show that 
V- 1 = V\ which implies [U\] a = V?* Q a p [U 3 k ]}. Ap- 
plying this relation to a Wilson line between -k (1> and 
-fc (2) , which are infinitesimally apart: 



viU QU [ulwb v^; Q aS (ii) 

=> W. fe(1 )^. fe(2) = V kW W* (1) ^ fe(2) V^ 2) (12) 



By employing (Bl|, the relation il2h is generalizable to 



finite-length paths, with arbitrary fc^ 1 , k (2> . We recall the 
definitions of C ky , T C ky , £ C ky and £TC ky , as detailed 



n App. [H] 

W\C ky ) = W{TC ky ) = V n , ky W{£TC ky y v_l tky 
= V^ ky W(C ky )*vl ky . (13) 



The first equality follows from ( B2 1 , the second from ( E3 ) 
and the third from £ T C ky = C~ Ky . Since V^ tk is unitary, 
we have shown that W along C ky is equivalent to the 
transpose of VV along C~ ky by a unitary transformation, 
thus proving (128]) . (|I3l) may be written in a familiar form: 



W -l ^-ky J = 0-1 W(£*» ) 6 (14) 

with = KVl k , and satisfying Q 2 = -I. This implies 
that each eigenstate of Wk at K y = {0, ir} has a degen- 
erate Kramer's partner. 

Appendix J: Isotropy of W 

Let us define {ip(k x )} ({$(k y )}) as Wannier trajecto- 
ries of the Wilson loop at constant k x (k y ). In this Section 
we prove that if {^(ky)} has a relative winding of W, so 
will {(p(k x )}. This follows because: 

(i) The Chern number is an obstruction to a smooth 
gauge in the BZ, hence both sets of trajectories, {^(ky)} 
and {ip(k x )}, must exhibit the same center-of-mass wind- 
ing Ci . 

(ii) Let M m (K x ), M W (K X ) and M (cc) (K x ) respectively 
be the number of +1,-1 and complex-conjugate-pair 
eigenvalues of the Wilson loop at constant K x = {0,7r}; 
M m (K x ) + M ( _ 1} (K X ) + M (CC) (K X ) = n occ , the number of 
occupied bands. We define 

M d = max{ Mf.^Tr) - M (4) (0) - M (cc) (0), 

M^W-M^W-M^W }; (Jl) 

if Mrf > there are number of Wannier trajectories 
that directly connect the atomic site (at K x ) and the 
mid-bond site (at K x +7r); if Md < there are none. By 
applying the mapping of Sec. |IIB[ it is possible to prove 
by exhaustion that Md — Nd, where Nd is defined in Eq. 
§. 

Given (i) and (ii) , the relations between W, C\ and Nd 
in Tab. |ll| imply that both sets of trajectories, {'ff(ky)} 
and {(^(fcT)}, have the same relative winding W. 
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